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Abstract 

We discuss microscopic mechanisms of the violent dynamics following strong laser excitation of a metal cluster 
embedded in a rare gas matrix, taking as test case Nag@Ar434. This covers at least two aspects: first, it represents 
the typical experimental situation of metal clusters handled in raregas matrices or a finite drop of surrounding 
raregas material, and second, it serves as a generic test case for highly excited chromophores in inert surroundings 
addressing questions of energy transport and perturbation of the medium. We simulate the process up to 10 ps 
using a mixed quantum mechanical (for the electrons) and classical (ions and atoms) approach and analyze the 
emerging dynamics with respect to all basic constituents: cluster electrons, cluster ions, and matrix atoms. We 
find several stages of relaxation taking place with time scales from a few fs to over a few ps and much slower 
processes remaining for long after the simulation. A particularly interesting aspect is that the surrounding raregas 
material stabilizes a highly charged metal cluster which would otherwise explode without delay. 
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1. Introduction 

Nanotechnology, which was predicted long ago 
to stimulate an industrial revolution [1] , now con- 
stitutes a fast developing research field, for exam- 
ple to increase the control of material at nanometer 
scale [2]. Especially metal clusters have attracted 
much attention in the past decades [3,4,5,6,7,8]. 
Among their many interesting characteristics and 
applications, their optical properties are particu- 
larly appealing [3]. These are dominated by the 
Mie surface plasmon which plays a key role in all 
dynamical regimes [8]. It provides a strong cou- 



pling to photons in a very narrow frequency win- 
dow which makes clusters an ideal laboratory for 
laser induced non-linear dynamics [9,10] . The plas- 
mon also serves as a versatile handle for pump and 
probe experiments [11,12,13,14,15], and it provides 
the key mechanism in driving the hefty Coulomb 
explosion of clusters [16,17]. For a recent summary 
on this great variety of phenomena of the dynam- 
ics of free metal clusters, see [18,19]. An even richer 
variety of phenomena emerges for metal clusters in 
contact with other materials, either embedded in 
a matrix or deposited on a surface. There is, e.g., 
second harmonic generation [20,21,22], dedicated 
shaping of clusters with intense laser pulses [23,24] , 
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and many more new phenomena emerging from the 
combination, for a review see [25]. Thus the proper 
combination of materials is a key issue in nano- 
technology [26]. There are also practical reasons 
for dealing with such compounds. The substrate 
simplifies the handling of the clusters such that 
many interesting experiments can only be done 
with clusters in contact with a carrier material 
[27,28,29,30]. Simple compoimds are also useful 
model systems for principle studies. For example, 
chisters in contact with insulators can serve as a 
versatile model system for chromophorcs which can 
be used as indicators in biological tissues [31,32] or 
for studies of radiation damage in materials [33,34] . 
The latter application deals with high energy den- 
sity related to intriguing questions of energy, mo- 
mentum and particle transport in materials. It is 
this aspect which we will address in the present 
manuscript. We are going here to study from a the- 
oretical perspective the (hindered) Coulomb ex- 
plosion of a metal cluster embedded in a raregas 
matrix. That as such is a topic of actual interest 
in cluster physics, see e.g [10,13,35,36,37]. And it 
serves as a well manageable test case for the more 
general issue of strong energy impact on matter 
and subsequent transport processes which are also 
of high current interest [38,39,40]. 

The theoretical description of embedded clus- 
ters is a great challenge because of the inevitable 
complications by the cluster-matrix interface and 
because of the huge number of atoms in the ma- 
trix. This holds even more so for the simulation 
of their dynamics. The quest for robust and ef- 
ficient schemes is omnipresent in all areas where 
one deals with large compounds. The solutions 
rely usually on a hierarchical modeling combining 
carefully quantum-mechanical and classical de- 
scriptions, see e.g. [41] for bio-molecules or [42,43] 
for chemical reactions on surfaces. Along these 
lines, we have developed a hierarchical model for 
the dynamics of metal clusters in contact with rare 
gas materials [44,45]. First successful tests of the 
performance of this approach in the moderately 
non-linear regime were presented in [46,47]. 

In this paper, we continue the study of the dy- 
namics of embedded metal clusters, considering 
now a violent laser excitation. The intensity is cho- 
sen such that a free cluster would explode imme- 



diately while the embedded configuration remains 

intact for a sufficiently long observation time of 
several ps. As test case, we consider the dynamics 
of a Nas cluster embedded in Ar434, following irra- 
diation by an intense and short (33 fs) laser pulse. 
We will pay particular attention to the impact 
of the environment on the dynamics of the metal 
cluster and analyze the interlaced processes from 
several points of view: detailed ionic and atomic 
trajectories, evolution of global shape parameters 
(radius, quadrupole deformation) with their asso- 
ciated vibrational time scales, and energy trans- 
fer providing information on the various relaxation 
times. Finally, we discuss a possible identification 
by pump and probe analysis in relation to recent 
experiments [23]. 



2. Formal framework 

In order to allow for a sufficiently large Ar ma- 
trix, we use a hierarchical approach. The many dif- 
ferent ingredients make a complete representation 
rather bulky. In that section, we sketch only briefiy 
the basic constituents of the approach and refer to 
[44,48] for a detailed layout. 

The Na cluster is treated in full micro- 
scopic detail using quantum-mechanical single- 
particle wavefunctions {(/?„(r, t),n = l...A^ei} for 
the valence electrons. These are coupled non- 
adiabatically to classical molecular dynamics 
(MD) for the Na ions which are described by their 
positions {R/,/ = l...A^Na}- The electronic wave- 
functions are propagated by time-dependent local- 
density approximation (TDLDA). The electron- 
ion interaction in the cluster is described by soft, 
local pseudo-potentials. This TDLDA-MD has 
been validated for linear and non-linear dynamics 
of free metal clusters [8,9]. 

Two classical degrees-of-freedom are associ- 
ated with each Ar atom: center-of-mass {Ro, a = 
1...Nai}, and electrical dipole moment which is 
parameterized as {R^ = Ro -|- do, a = l...iVAr}- 
Note that the Ar dipole is practically handled by 
two constituents with opposite charge, positive Ar 
core (at Ro) and negative Ar valence cloud (at 
Ro). With the atomic dipoles, we explicitely treat 
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the dynamical polarizability of the atoms through 
polarization potentials [49]. Smooth, Gaussian 
charge distributions are used for Ar ionic cores 
and electron clouds in order to regularize the sin- 
gular dipole interaction. We associate a Gaussian 
charge charge distribution to both constituents 
having a width of the order of the 3p shell in Ar, 
similar as was done in [51]. The Coulomb field of 
the (softened) Ar dipoles provides the polarization 
potentials which are the dominant agents at long 
ranges. The Na"*" ions of the metal cluster have 
net charge ^Na = +1 and interact with the Ar 
dipoles predominantly by the monopolc moment. 
The small dipole polarizability of the Na+ core 
is neglected. The cluster electrons do also couple 
naturally to the Coulomb field generated by the 
atoms. 

The model is fully specified by giving the total 
energy of the system. It is composed as 

-E'total = -E^Nacluster + ^^Ar + ^^coupl + -E'vdW • 

The energy of the Na cluster -BNaciustor consists out 
of TDLDA for the electrons, MD for ions, and a 
coupling of both by soft, local pseudo-potentials; 
that standard treatment is well documented at 
many places, e.g. [8,9]. The Ar system and its cou- 
pling to the clusters are described by 
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The calibration of the various contributions is 
taken from independent sources, except perhaps 
for a final fine tuning to the NaAr dimer modify- 
ing only the term WeiAr of Eq.(4). The parameters 
are summarized in table 1. The third column of 
the table indicates the source for the parameters. 
In the following, we report briefly the motivations 
for the choices. 

Most important are the polarization potentials 
defined in Eq.(3). They are described by the model 
of a valence electron cloud oscillating against the 
raregas core ion. Its parameters are: gAr the effec- 
tive charge of valence cloud, rriAr = (/ArWei the 
effective mass of valence cloud, fcAr the restoring 
force for dipoles, and ctaf the width of the core and 
valence clouds. The gAr and fcAr are adjusted to 
reproduce the dynamical polarizability aoiu)) of 
the Ar atom at low frequencies, i.e. we choose to 
reproduce the static limit an (w = 0) and the sec- 
ond derivative of a'lj{J' =0). The width dAr is de- 
termined consistently such that the restoring force 
from the folded Coulomb force (for small displace- 
ments) reproduces the spring constant fcAr- 

The short range repulsion is provided by the var- 
ious core potentials. For the Ar-Ar core interaction, 
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Table 1 

Parameters for the various model potentials. 

Eq.(5), we employ a Lennard- Jones type potential 
with parameters such that binding properties of 
bulk Ar are reproduced. The Na-Ar core potential, 
Eq.(6), is chosen according to [50]. Note that the 
Na-Ar potential from [50] does contain also a long 
range part oc uat which accounts for the dipole 
polarization-potential. We describe that long range 
part explicitely and have to omit it here. We thus 
choose the form as given in V^rNa- 

The pseudo-potential WeiAr, Eq.(4), for the 
electron-Ar core repulsion has been modeled ac- 
cording to the proposal of [51]. Its parameters 
determine sensitively the binding properties of Na 
to the Ar atoms. We use them as a means for a 
final fine-tuning of the model. The benchmark for 
adjustment is provided by the Na-Ar dimcr. The 
data (dimer binding energy, bond length, exci- 
tation energy) are taken from [60] and [61]. The 
adjustmcnct shows some freedom in the choice 
of Aei. We exploit that to produce the softest 
reasonable core potential. 

The Van-der-Waals energy ii^vdw, Eq.(8), is a 
correlation from the dipole excitation in the Ar 
atom coupled with a dipole excitation in the clus- 
ter. We exploit that lomic is far below the excita- 
tions in the Ar atom. This simplifies the term to 
the variance of the dipole operator in the cluster, 
using again the regularized dipole operator de- 
fined in Eq. (9), corresponding to the smoothened 
Ar charge distributions. The full dipole variance is 
simplified in terms of the local variance. 

The classical MD for Na ions, Ar atoms and Ar 
dipoles is treated in full three dimensions. How- 
ever, the quantum mechanical treatment of the Na 



cluster involvc;s two approximations. The Kohn- 
Sham equations for the electron cloud involves ax- 
ially averaged potentials and the self-interaction 
correction (SIC) to TDLDA is treated at the level 
of average-density SIC (ADSIC) [62,63]. Both ap- 
proximations enhance somewhat the barriers for 
fragmentation of the cluster. The axial approxima- 
tion and ADSIC are needed to allow the large scale 
calculations performed here. The SIC, on the other 
hand, is compulsory for an appropriate description 
of ionization which is the leading process in the ini- 
tial stage of laser excitation. These two approxima- 
tions mean that we provide rather a qualitative pic- 
ture. The effects shown are certainly relevant. The 
thresholds concerning charge state and laser field 
strength are probably not too precise. However we 
also dispose of a full three-dimensional treatment 
for the valence electrons which, moreover, can han- 
dle various levels of SIC and even exact exchange. 
We used that for occasional counterchecks and we 
never found significant deviations from the approx- 
imate treatment. 

Finally, a few words are in order about the pos- 
sible range of validity for that hierarchical model. 
The structure of the coupled systems is well de- 
scribed by construction, see the predecessor in [51] 
and the extensive testing in [44] . The same holds for 
optical response (sec [45,64]). However, one has to 
remain aware that the model has limitations with 
respect to allowed frequencies and amplitudes. The 
dynamical response of the Ar atoms is tuned to fre- 
quencies safely below the Ar resonance peak, i.e. 
safely below 1 Ry. This is well fulfilled in all our 
calculations because the leading frequency is the 
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cluster's Mie plasmon at around 0.2 Ry. The am- 
plitudes of the dipole oscillations in Ar should stay 
below the threshold where free electrons are re- 
leased from the Ar atom into the matrix. We have 
checked that process by fully quantum mechanical 
calculations of laser excitation on Ar atoms and we 
find a critical field strength of order of 0.1 Ry/ao. 
We protocol during our calculation the actual field 
strength at each Ar site and thus check that we do 
not exceed that limit in our simulations. 



3. The test Ccise 

We focus here on a Nag cluster embedded in a 
large Ar cluster with 434 atoms. This should serve 
BjS Bj finite model of an Ar matrix, although as- 
pects of finite mixed clusters will also be addressed 
briefly. From a more general perspective, we en- 
counter here a typical example of a chromophore 
inside an inert environment. The case is particu- 
larly interesting because the mix of simple mate- 
rials allows to concentrate on the basic effects of a 
chromophore without the need to bother about the 
elaborate technical details of more complex mate- 
rials, as e.g. large organic molecules. We consider 
the case of low initial temperatures, safely below 
the Ar melting point of 84 K, where spatial fluctua- 
tions of the atoms and ions can be safely neglected. 

The small Na cluster at low temperature has a 
very clean excitation spectrum with rather narrow 
excitation peaks in the plasmon range around 2.5 
eV [8] . Correlation eifects beyond TDLDA remain 
negligible for the global observables considered 
here. Spectral fragmentation from coupling to 
energetically close electron-hole states (Landau 
fragmentation) is negligible for that system size 
because the Mie surface plasmon resides in a gap 
of such states [53] . Broadening from higher corre- 
lations is also small [54] , again due to the spectral 
gap for that magic electron number. Thermal 
broadening in the initial state is negligible at the 
presently studied low temperatures [52]. Larger 
electronic temperatures develop in the course of 
the dynamics and they call, in principle, for a 
description beyond pure mean field. Electron- 
electron collisions could be included by switching 



to a semi-classical Vlasov-Uehling-Uhlenbeck de- 
scription of the electron cloud [48,65]. But long- 
time semi-classical propagation is still extremely 
hard to stabilize and, as we will see in section 4.3, 
electrons are thermally almost decoupled from 
ions. Thus we can neglect that effect for time scales 
explored in the present study. The surrounding Ar 
cluster has to be sufficiently large to absorb the 
energy from the highly excited chromophore with- 
out breakup and to avoid specific size or surface 
effects. We have checked matrix sizes from 30 to 
f 048 Ar atoms for the dynamical regime intended 
here and found qualitatively similar behaviors 
from f 64 Ar atoms on. To be on the safe side, the 
present results have been obtained for a model 
matrix of 434 Ar atoms. The setup was generated 
by cutting a closed pack of 447 Ar atoms out of 
a crystal, removing further 13 atoms from the 
center, placing the Nas in the cavity, and finally 
re-optimizing the whole configuration by Monte- 
Carlo simulated annealing techniques. The Nag 
cluster with a radius of 7 ao (r.m.s. -radius of 5.5 
ao) then resides in a cavity of about 10 ao radius. 
The outermost Ar shell has a radius of 31 ap. This 
provides a substantial coverage of the embedded 
Nag. The spatial structure is still very close to 
the configuration before re-optimization [45] . The 
Nag consists out of two rings of four ions each and 
the rings have a relative rotation of 45° to mini- 
mize the Coulomb repulsion between the ions. The 
global shape is axially symmetric, only slightly 
oblate deformed. The surrounding Ar maintains 
basically its original crystal structure with some 
small rearrangement for the innermost shells. 

The dynamics in the linear regime was studied 
in [45]. A first case of moderate laser excitations 
was discussed in [46] . Here we aim at a more vio- 
lent, still non-destructive, process. To that end, we 
irradiate the Na8@Ar434 system with a laser of fre- 
quency 1.9 eV, intensity 2.4 * 10^^ W/cm^, and a 
cos^ pulse envelope with total pulse length of 100 
fs (FWHM = 33 fs). The laser is polarized along 
the symmetry axis of the system, henceforth called 
z axis. The rather short and intense laser excita- 
tion produces very quickly a net ionization stage 
3+ and deposits a large amount of energy in the 
system. We then continue to propagate the system 
up to about 10 ps in order to cover the relevant 
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time scales. 



4. Results and discussion 

4.1. Detailed view at the propagation of ions and 
atoms 

The initial reaction of the system is dominated 
by the large dipole response of the electrons. For 
the short pulses used here, the strong electrical field 
with subsequently large dipole amplitude leads to 
a direct emission of 3 electrons which escape in- 
stantaneously (i.e. within 3-10 fs). The thus pro- 
duced large Coulomb pressure leads to a rearrange- 
ment of the whole system (Na ions and Ar atoms) 
at ionic/atomic time scale as is demonstrated in 
figure 1 showing the time evolution of Na ions and 
Ar atoms in detail. It exhibits different int(T(^st- 
ing time scales. The first time scale (after the elec- 
tronic one, of order of a few fs) is observed in the 
reaction of the Na ions. They perform a Coulomb 
explosion whose first signs can be seen at about 
50 fs. That proceeds almost identical to the sim- 
ilar case of a free Nag cluster lifted to charge 3+. 
But the "explosion" is abruptly stopped at about 
200 fs when the ions hit the repulsive core of the 
first shell of Ar. The ionic motion turns to damped 
oscillations around a (r.m.s.) radius of about 7 ao. 
In atoms, one also observes two stages, although 
on longer time scales. The first phase is a spread- 
ing of the perturbation (caused by the stopping of 
Na ions) into the various Ar shells. This is espe- 
cially visible along the z axis which allows to read 
off the propagation speed of this perturbation as 
20-30 ao/ps. We have estimated the sound velocity 
in the corresponding large pure Ar cluster (Ar447) 
by computing its vibrational spectrum. The radial 
compression mode corresponds to the sound mode 
in bulk material. Its frequency is Wyib = 1-8 meV. 
The momentum of the radial wave is q = w/R 
where i? = 30 ao is the cluster radius. The sound 
velocity is then estimated as t^sound = ^vih/q ~ 30 
ao/ps which is very close to the propagation speed 
as observed in the figure. This suggests an interpre- 
tation as a sound wave sent by the initial bounce 
of the Na ions. 



Once transferred to a given Ar shell the pertur- 
bation generates oscillations combined with some 
diffusion which after about 1.5 ps has spread over 
all shells. The perturbations affect even the outer- 
most shell which thus acquires some oscillations. 
That is an effect of finite size of the "matrix" . The 
oscillation amplitude shrinks with increasing sys- 
tem size and grows towards smaller matrices. For 
example, it was found to be smaller for the larger 
test case with 1048 Ar atoms and the last Ar atoms 
are even kicked away for matrices at and bellow 164 
Ar atoms. The reason is twofold: smaller system 
provide, first, less Ar shells to absorb energy, and 
second, less dipole binding as whole. This question 
will be discussed more systematically in a larger 
forthcoming paper. 

The relaxation of matrix oscillations is much 
slower than that of Na ions and beyond the time 
scale computed here. These long time scales for full 
relaxation and evaporation of Ar atoms are well 
known from experiments of dimer molecules em- 
bedded in Ar clusters, see e.g. [55] . 

We finally comment briefly on the flnite net 
charge of our test system. It emerges because the 
electrons propagate almost unhindered through 
the surrounding Ar cluster and eventually escape 
to infinity. In a macroscopically large matrix, the 
electrons would be stuck somewhere within the 
range of their mean free path, drift very slowly 
back towards the now attractive cluster well, 
and eventually recombine there. We have simu- 
lated that by using reflecting boundary conditions 
rather than absorbing ones and we find recombi- 
nation times of the order of several 10 ps. Further- 
more the ionic/atomic rearrangements following 
the irradiation turn out to be qualitatively very 
similar whatever the boundary conditions. Thus 
the present scenario should provide a pertinent 
picture for several ps, the time window studied 
here. 

4.2. Time evolution of global shapes 

Figure 2 visualizes the ionic/atomic dynamics 
in terms of simple global observables. The middle 
panel shows the ionic r.m.s. radius which provides 
an even clearer picture of the ionic dynamics de- 
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Fig. 1. Time evolution of the the atomic (full lines) and ionic (dotted lines) z-coordinates (lower panel) and radial distances 

r = y^r^ + j/^4-z^ (upper panel) for Na8Ar434 excited with a laser of frequency = 1.9 eV, intensity 2.4 * 10^^ W/cm'^, 
and a cos"^ pulse profile with FWHM = 33 fs. The laser was polarized along the z axis which is also the symmetry axis of 
the system. 



scribed above. Let us first concentrate on the the 
lower two panels showing the global radial oscil- 
lations (breathing). The electronic radius expands 
during the short excitation phase. However, as soon 
as the laser is switched off after 100 fs, the elec- 
tronic radius quickly relaxes towards the ionic one, 
as can be seen in the lower panel of figure 2. Indeed 
the strong monopole Coulomb forces ties electronic 
and ionic radius always closely together. The ra- 
dial oscillations of the cluster have the typical cy- 
cle of about 250 fs, well known from free Na clus- 



ters [56] and re-established for embedded ones [46] . 
The most interesting aspect is here the clearly vis- 
ible relaxation of the cluster oscillations at a time 
scale of about 4 ps. It seems that a slow energy 
transfer from the ionic cluster vibrations into the 
raregas environment takes place. This will be in- 
spected further in section 4.3. The breathing os- 
cillation of the Ar environment (i.e. of the atomic 
radii) proceeds at a cycle of about 3 ps (see mid- 
dle panel of figure 2), corresponding to a volume 
frequency of about 1.4 meV. This is within fac- 
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Fig. 2. Time evolution of global shape observables of atoms (solid lines), ions (dotted lines), and electrons (fine dashed 
lines). Lowest panel: r.m.s. radii of cluster constituents, ions and electrons. Middle panel: r.m.s. radius of Ar atoms. Upper 
panel: Quadrupole moments of ions and atoms. 



tor two the surface vibration frequency observed in 
neutral Ar clusters of that size [57] , which is quite 
gratifying in view of the structural differences (Na 
cluster embedded, high net charge) between these 
two systems. A detailed comparison will have to 
be worked out in future studies. 

The upper panel complements the analysis by 
the quadrupole moments of the Na cluster (ions) 
and of the Ar atoms. The quadrupole oscillations 
of both subsystems arc again at their very own 
time scale which a very long cycle time, about 7 
ps, for quadrupole motion of the Ar atoms. The Na 
cluster shows, moreover, a global trend to finally 



oblate shapes. Note that a similar effect, the build- 
up of persistent quadrupole deformations following 
laser irradiation, has been observed for Ag clus- 
ters embedded in glass [23] . Glass is a much harder 
material than Ar but we may speculate that the 
strong laser pulse provides also an intense heating, 
and thus melting, of the immediate neighborhood 
of the embedded cluster. For then the mechanisms 
of deformation may be similar. 
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Fig. 3. Lower left; time evolution of kinetic energies for atoms and ions; the line denoted "non-coUcctivc" shows the internal 
kinetic energy of the atoms after subtracting the contribution from global radial breathing. Upper left: Distribution of atomic 
kinetic energies at time t = 9.5 ps where the boxes show the actual atomic distribution sampled in energy bins and the dotted 
line stands for the equivalent thermal (MaxwcUian) distribution with the temperature as given by the instantaneous kinetic 
energy. Right: energy balance, i.e. input through the laser field (left) and its distribution over various subsystems (right). 



4.3. Energetic aspects 

The lower left panel of figure 3 shows the kinetic 
energies of ions and atoms. The Na ions gather en- 
ergy very quickly from the starting Coulomb ex- 
plosion. At around 200 fs, the quickly acquired en- 
ergy is transferred abruptly to the atoms, leaving 
only one fifth of kinetic energy for the first phase of 
ionic oscillations. Prom then on, both kinetic ener- 
gies proceed for a while at nearly constant values 
with some residual oscillations. At around 2 ps, a 
slow transition seems to take place where the ionic 
energy flows towards the atomic subsystem. The 
situation stabilizes after 3 ps to a thermal share of 
energies where practically all kinetic energy resides 
in the atoms and only a tiny fraction remains for 
the ions. However, the notion of "thermal share" 
does not imply a full equilibrium even at final sim- 



ulation time of 9.5 ps. We know from figure 2 that 
there are very slow and long lasting atomic vibra- 
tions. We should call the present situation more 
correctly a local (or instantaneous) equilibrium. 

The global atomic breathing oscillations carry 
still a large amount of energy. To quantify that 
aspect we have computed the collective energy. The 
fine dashed line in the lower left panel shows the 
remaining "internal" kinetic energy of the atoms. 
And indeed, almost two thirds of the energy reside 
still in the collective breathing yet awaiting final 
thcrmalization at a time scale far beyond 10 ps. 
Nonetheless, it is interesting to have a look at the 
distribution of kinetic energies. This is shown in the 
upper loft panel of figure 3. It is the distribution of 
the total kinetic energies and it is compared with 
the thermal distribution at a temperature of 20 
K, equivalent to the full kinetic energy at the late 
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times in the plot. Both distributions compare fairly 
well in view of the fact that full equilibrium has 
not yet been reached (thus allowing for somewhat 
larger fluctuations). 

The right panel sketches the main energy share 
in the reaction. The system acquires as much as 23 
eV energy from the laser pulse. The by far largest 
portion, namely 83%, is carried away very quickly 
by the emitted electrons. Two third from what re- 
mains flow into ionic and atomic motion. We have 
learned from the lower left panel that this encTgy 
arises from the stopped Coulomb explosion, and 
is shared in two steps between ions and atoms: 
initially, the ions acquire a large potential energy 
from sudden charging, this is converted through 
"Coulomb explosion" instantaneously to a large 
kinetic energy, which then is transferred quickly 
to the atoms when the "explosion" is abruptly 
stopped. A somewhat slower second phase of en- 
ergy transfer follows after 2 ps which is related to 
rearrangement and some exchange with still avail- 
able potential energy. One third of the remain- 
ing energy, i.e. about 1 eV, remains with the elec- 
trons as internal kinetic energy (this is the excess 
kinetic energy beyond the unavoidable quantum 
mechanical minimum, see [9]). That is very little 
in view of the fact that it has been the cluster's 
electron cloud which served as doorway to absorb 
all energy from the laser field. It is, however, a 
great deal in terms of temperature. The electronic 
temperature corresponds to 1500 K whereas the 
ionic/atomic temperature has yet come up to only 
18 K. That makes it clear that the process resides 
still in a transient stage where the electronic and 
the ionic/atomic subsystems have reached a local 
equilibrium amongst themselves. Full equilibration 
and subseriucnt evaporative cooling will follow at a 
much longer time scale, far beyond the capabilities 
of the very detailed description through a TDLDA- 
MD approach. 

There remains the question which energetic as- 
pects are generic for metal clusters in a rare gas en- 
vironment and which are specific to the particular 
mixed system studied here. The generic aspects are 
related to what happens with the cluster and its 
immediate rare gas vicinity. This concerns the en- 
ergy absorption from the laser which proceeds pre- 
dominantly through the metal cluster's electrons, 



the stopping of the Coulomb explosion and associ- 
ated fast energy transfer to the atoms, the relax- 
ation time scales in these early phases of a few ps, 
and the energy balance in absolute numbers. The 
particular property which depends on the matrix 
size is most of all the atomic temperature. We find 
that the same absolute amount of energy is trans- 
ferred from the laser to the Nag cluster, indepen- 
dent of matrix size. The energy loss through ion- 
ization is also the same. Thus the final atomic tem- 
perature will then be inversely proportional to the 
matrix size. 

4.4. Possible pathways to experimental analysis 

The question remains how one could analyze the 
key pattern of the embedded cluster dynamics ex- 
perimentally. The most interesting effect is the hin- 
dered explosion of the imprisoned Na cluster and 
its subsequent shape oscillations in the Ar cavity. 
Such shape oscillations and eventual relaxation to 
deformed shape have been produced and observed 
experimentally for Ag clusters embedded in glass 
[23] or deposited on a substrate [58]. The (time- 
dependent) cluster deformation is verified by probe 
pulses measuring the actual optical response of the 
metal cluster, exploiting the pronounced features 
of the Mie plasmon resonance. Indeed, the plasmon 
couples strongly to light at very specific frequen- 
cies and these plasmon frequencies are uniquely 
related to the cluster shape. For free clusters, we 
have discussed a rather simple pump-and-probe 
setup which allows to map in a unique fashion the 
evolution of radial shape [12], of quadrupole de- 
formations [59] and of elongation in a fission pro- 
cess [14,15]. We have computed that for the case 
considered here. To this end, we take the instanta- 
neous configuration at a given time and compute 
the optical response in all three spatial directions. 
These expensive calculations have been performed 
in full 3D TDLDA (which allows to properly ac- 
cess the directions perpendicular to the laser polar- 
ization axis) on a representative time of extremal 
deformations. The other points are interpolated 
by well benchmarkcd plasmon calculations for a 
smaller subsystem. Results are shown in figure 4. 
The radial oscillations (middle panel of figure 2) re- 
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Fig. 4. Time evolution of the estimated Mie plasmon fre- 
quencies in rr, y and z directions for the embedded Nas 
cluster. The cluster remains always almost axially symmet- 
ric such that the x and y modes are nearly degenerate. The 
spectral distribution at the last time point is indicated to 
the right of the figure. 

late to the oscillations of the average peak position 
which are particularly well visible here in the early 
phase up to about 1.5 ps. The then evolving oblate 
deformation leads to a splitting of the resonance 
peak where the shorter extension along z is associ- 
ated to a blue shift of the mode and the larger ex- 
tension in orthogonal direction to a red shift. The x 
and y modes are nearly degenerate due to near ax- 
ial symmetry of the total system. Pump and probe 
analysis of the time-dependent plasmon peaks thus 
provides worthwhile information about the shapes 
of the embedded cluster, which can be measured 
by probe pulse as explained in [12,59,14,15]. Note 
here that the experimental results obtained for Ag 
clusters in glass [23] are in qualitative agreement 
with our findings. Our model thus provides a mi- 
croscopic interpretation, the first one, to the best 
of our knowledge, of these experimental investiga- 
tions. 



5. Conclusion 

In this paper, we have analyzed the dynamics 
of a Na cluster embedded in an Ar matrix follow- 
ing irradiation by an intense laser pulse. We have 
used a hierarchical model with a detailed quantum- 
mechanical treatment of the electron cloud in the 
Na cluster and a classical description of the Na 



ions and surrounding Ar atoms whereby we include 
the Ar dipole moments to account for the dynam- 
ical polarizabilities. As test case, we have consid- 
ered Na8@Ar434 which provides several interest- 
ing aspects : It is a prototype of a chromophore 
in an else-wise inert material, it serves as model 
for a cluster embedded in a raregas matrix, and 
wc do also find interesting features which arc spe- 
cific for a finite mixed cluster. For the excitation 
mechanism, we consider irradiation by a short 33 
fs laser pulse with intensity of about lO^^W/cm^. 
The laser couples almost exclusively to the clus- 
ters electron cloud and deposits about 23 eV en- 
ergy in the system. This leads to immediate ioniza- 
tion of Nag to a 3+ charge state. The subsequent 
Coulomb explosion of the Nag is abruptly stopped 
by the Ar atoms, turning the motion into steady 
shape oscillations of the cluster while the absorbed 
momentum spreads in a radial sound wave prop- 
agating through the atoms and triggers global ra- 
dial oscillations thereof. A second phase around 2- 
3 ps transfers most of the remaining ionic energy 
to the atoms and cools the ions to almost atomic 
temperature. The relaxation process can be nicely 
seen in oscillations of cluster radius and deforma- 
tion. Both these observables are accessible to ex- 
perimental observation by an especially designed 
pmnp-and-probe setup. 

A further interesting feature is that the Ar sur- 
roundings do not only stop the Coulomb explosion. 
They even stabilize the high charge state at a time 
scale which is sufficiently long for an observation 
of ionic and atomic motion. This is probably a par- 
ticular feature of finite mixed clusters. The sur- 
rounding Ar system should be large enough to de- 
liver stabilization through induced dipole attrac- 
tion, but small enough such that the cluster elec- 
trons which are directly emitted due to the laser 
excitation can travel nearly unhindered through 
the Ar shells and finally escape to whole system. 
The degree of charge stabilization in dependence 
of system size and the mean free path of electrons 
through ordered raregas shells is a task for future 
studies. 
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